 Lector univ dr Cristina Nartea Cursul 4 Matrice Rangul unei matrice Rezolvarea sistemelor de ecuatii liniare Metoda eliminarii a lui Gauss Definitie O matrice m × n este o serie de mn intrari, numite elemente, aranjate in m linii si n coloane in cazul in care o matrice se noteaza cu A, elementul din randul i si coloana j se noteaza cu aij si matricea se scrie  a11  a1n    A       a   m1  amn  Exemple de matrice  Matrice 1x3 1 2 3  2  Matrice 3x1   5 6    2 4  Matrice 2x2   1 0  1 0 2    Matrice 3x3  1 3 4   5 1 2   O matrice patratica este o matrice in care numarul de linii m este egal cu numarul de coloane n Egalitatea a doua matrice Egalitatea a doua matrice inseamna ca, daca A si B sunt egale, atunci fiecare este o copie identica a celeilalte 1 2 3 1 2 3 Ex A    si B    Aflati x astfel incat A=B  0 1 2  0 x 2 Adunarea a doua matrice Adunarea de matrice A si B este definita numai in cazul matricele au acelasi numar de randuri si cu acelasi numar de coloane Sa consideram A  a i j  si Lector univ dr Cristina Nartea B  bi j  sa fie matrice m × n Matricea m × n formata incat elementul din linia i si coloana j este aij  bij pentru fiecare i si j este matricea A + B 1 2  1 0  2 2 Ex Pentru A   , B   , A B    3 4  1 1  2 5 inmultirea unei matrice cu scalari Sa consideram matricea m × n, A  a i j  si λ un scalar (real sau complex) in cazul in care A este inmultit cu λ, si se scrie λA, fiecare element din A este inmultit cu λ pentru a obtine matrice m × n,  A  a i j  1 2 3 2 4 6 Ex Pentru  =2 si A   , A     0 1 2   0 2 4  inmultirea matricelor Este important sa observati ca, atunci cand produsul AB este definit, produsul BA este in general diferit sau poate sa nici nu fie definit Se pot inmulti matrice de tip mxn cu matrice nxp, iar rezultatul este o matrice de tip mxp 1  1 0 2   Ex A    , B   2  Rezultatul este  1 4 5   3    1 1  0  2  2  3   7  A B        (1) 1  4  2  5  3   22   1 5  1 0 1 0 1   Tema 1 Fie A   , B    , C    , D   4 2  Determinati care inmultiri pot fi  2 3   0 1     3 1 0    efectuate si in acest caz calculati: AB, BA, AC, CA, ABC, CAB, AD, DA, CD, DC, ACD, DAC Transpusa unei matrice Sa consideram matricea m × n, A  a i j  Atunci transpusa lui A, notata de A T este matricea obtinuta schimband liniile in coloane pentru a produce o matrice nxm, AT   a ji  1 2 T  1 3 Ex A   , A    3 4  2 4 Lector univ dr Cristina Nartea  1 2  1 0 1 T   A , A   0 1 2 1 3   1 3    Determinantul unei matrice Fiecare matrice patratica, are ca element asociat un singur numar determinant al lui A Daca A este o n ×n matrice, determinantul lui A este indicat prin afisarea elementelor lui A intre doua bare verticale, dupa cum urmeaza: a11  a1n    an1  ann 1 2 Ex Determinant de ordinul 2  1 (4)  2  5  14 5 4 Determinant de ordinul 3 1 0 1 2 3 4  1 3  6  2  2 1  (1)  0  4  1 3  (1)  4  2 1  6  0  2  18  4  3  8  17 1 2 6 1 0 1 2 3 4 Rangul unei matrice Fie A  M m,n    o matrice nenula Spunem ca matricea A are rangul r si notam rang A  r , daca A are un minor nenul de ordin r, iar toti minorii lui A de ordin mai mare decat r (daca exista) sunt nuli Aplicatie Calculati rangul matricelor  1 0 2  1 0   B   1 3 C  1 0 2   A   1 3 0 ,  ,  , D   3 4  2 2 1  0 4 3 1 4  1 2      1 inversa unei matrice Daca det A  0 , atunci A este inversabila si A1   A* det A 1 2 Ex Calculati inversa matricei A    3 1 Lector univ dr Cristina Nartea 1 2  det A   1  6  5  0 , deci A este inversabila 3 1  1 3  AT     2 1  Complementii algebrici 11  (1)11 1  1 12  (1)1 2  2  2  21  (1)21  3  3  22  (1)2 2 1  1  1 2   A*     3 1   1 2   1 1  1 2   5 5   A     5  3 1   3 1     5 5 1 0 2 Aplicatie Calculati inversa matricei A   2 1 3  0 1 4   Aplicatie la spatii vectoriale Fie B={ e1  (1,1), e2  (2,3) } si B’={ f1  (1,3), f 2  (3,8) } a) Sa se verifice daca sistemele de vectori B si B’ formeaza baze b) Sa se gaseasca matricea de trecere din baza B in baza B’ c) Daca x[ B ']  (2,1) , gasiti coordonatele lui x in baza B d) Daca x[ B ]  (1, 0) , gasiti coordonatele lui x in baza B’ Sisteme de ecuatii liniare Forma generala a unui sistem de m ecuatii liniare cu n necunoscute este: a11 x1  a12 x2    a1n xn  b1 a x  a x    a x  b  21 1 22 2 2n n 2  (1 1)  am1 x1  am 2 x2    amn xn  bm unde:  x1 , x2 , xn sunt necunoscutele sistemului, Lector univ dr Cristina Nartea  numerele aij , i  1, m, j  1, n sunt coeficientii necunoscutelor,  b1 , b2 ,bm sunt termenii liberi ai sistemului Unui sistem liniar ii asociem urmatoarele matrice:  a11 a12  a1n    a a 22  a 2n   A   21 matricea sistemului,        a  a mn   m1 a m 2  b1     b2      matricea termenilor liberi   b   m  x1     x2      matricea necunoscutelor,   x   n  a11 a12  a1n b1       a 21 a 22  a 2n b2   A matricea extinsa a sistemului care se obtine        a  a mn bm   m1 a m 2 adaugand la matricea A coloana termenilor liberi Definitia 1 Se numeste solutie a sistemului de ecuatii liniare un sistem ordonat de n numere 1, 2 , n t astfel incat inlocuind necunoscutele x1 , x2 , xn respectiv prin 1 ,  2 ,  n este verificata fiecare din ecuatiile sistemului Definitia 2 Un sistem este  compatibil daca are cel putin o solutie,  compatibil determinat daca are solutie unica,  compatibil nedeterminat daca are o infinitate de solutii,  incompatibil daca nu are solutii Metode de rezolvare a sistemelor liniare Lector univ dr Cristina Nartea 1) Metoda lui Cramer permite rezolvarea sistemelor liniare de n ecuatii cu n necunoscute avand determinantul asociat matricei sistemului nenul Teorema 1 Daca sistemul a11 x1  a12 x2    a1n xn  b1 a x  a x    a x  b  21 1 22 2 2n n 2  (1 2)  an1 x1  an 2 x2    ann xn  bn are determinantul  nenul, atunci solutia sa utilizand metoda lui Cramer este  x1 ,, xn  , unde xi xi  , i  1, n , xi , i  1, n fiind determinantul obtinut din  prin inlocuirea coloanei  corespunzatoare coeficientilor necunoscutei xi , i  1, n cu coloana termenilor liberi, adica a11 a12  a1,i 1 b1 a1,i 1  a1n a 21 a 22  a 2,i 1 b2 a 2,i 1  a 2 n xi          a n1 a n 2  a n ,i 1 bn a n ,i 1  a nn 2) Metoda de rezolvare a sistemelor liniare de m ecuatii cu n necunoscute 1) Se determina rang A a11 a12  a1r a21 a22  a2 r 2) Se alege un minor principal  p      ar1 ar 2  arr 3) Se precizeaza: necunoscutele principale x1 ,, xr si secundare xr 1 , xr  2 , xn si de asemenea ecuatiile principale (ecuatiile 1, 2, r ) si ecuatiile secundare (celelalte m  r ecuatii) Daca exista ecuatii secundare se calculeaza minorii caracteristici (minorul obtinut din minorul principal, prin bordarea acestuia cu elementele corespunzatoare ale coloanei termenilor liberi si cate una din liniile ramase); numarul minorilor caracteristici este egal cu numarul ecuatiilor secundare si este egal cu m  r Lector univ dr Cristina Nartea 4) Se stabileste daca sistemul (1 1) este compatibil Teorema 2 (Teorema lui Rouche) Un sistem de ecuatii este compatibil daca si numai daca toti minorii caracteristici sunt nuli Teorema 3 (Teorema Kronecker - Capelli) Conditia necesara si suficienta ca sistemul sa fie compatibil este ca rangA  rangA 5) Daca sistemul este compatibil solutia sa se obtine prin rezolvarea sistemului principal format din ecuatiile rezultate trecand in membrul drept termenii care contin necunoscutele secundare si atribuind acestor necunoscute secundare valori arbitrare): - daca numarul necunoscutelor secundare este 0 sistemul este compatibil determinat; - daca exista necunoscute secundare, sistemul este compatibil nedeterminat; numarul necunoscutelor secundare arata gradul de nedeterminare 3) Metoda transformarilor elementare (Metoda eliminarii a lui Gauss) Metoda transformarilor elementare este de fapt procedeul de reducere a necunoscutelor, scris, eventual, sub forma matriceala in cazul sistemelor de doua ecuatii cu doua necunoscute, aceasta metoda este de fapt metoda reducerii Exista 3 tipuri de transformari elementare Schimbarea a doua ecuatii; inmultirea unei ecuatii cu un scalar nenul; Adunarea unei ecuatii inmultite cu un scalar la o alta ecuatie 2 x  3 y  4 Exemplul 1 Rezolvati sistemul   x  y  5 Sistemul Matricea extinsa si transformarile elementare 2 x  3 y  4  2 3 | 4     x  y  5  1 1 | 5  Lector univ dr Cristina Nartea  1 L  L 2 2 1 2 x  3 y  4 2 3 | 4    5  0 5  2 y  7 | 7   2   42  22  11 2 x  5  4  2x    x  5  5     y  14  y  14  y  14  5   5  5 x  3y  z  3  Exemplul 2 Rezolvati sistemul  x  y  2 z  2 3x  2 y  z  1  Solutie x  3y  z  3  1 3 1 | 3      x  y  2 z  2  1 1 2 | 2  3x  2 y  z  1  3 2 1 | 1     L2  L1  L3 3L1 x  3y  z  3  1 3 1 | 3       2 y  3z  5  0 2 3 | 5   11y  4 z  8  0 11 4 | 8      11 L3  L 2 2 Lector univ dr Cristina Nartea    x  3y  z  3  1 3 1 3  |      2 y  3z  5  0 2 3 | 5   25 39  25 39   z 0 0 |   2 2  2 2   39  24  x  3 y   3  x  25 25    39  4   2 y  3  5   y    25  25  39  39  z  z  25  25 Calculul inversei unei matrice prin metoda transformarilor elementare Aplicatie Sa se determine inversele matricelor  2 3 a) A     5 1  2 1 1  b) A   1 2 3   3 1 1    2 3 | 1 0  2 3 | 1 0  L2 L2  52 L1    2 L2   L2 Solutie a)      13 5 1 | 0 1  0  13 |  5 1  2 2   2 6   1 3  2 3 | 1 0   2 0 |    1 0 |  13 13  1      L1  L1 3 L2  13 13 L   1 2 L1    0 1 | 5 2   0 1 | 5 2  5 2  13 13     0 1 |    13 13   13 13   1 3    13 13  Deci A1     5 2     13 13  Lector univ dr Cristina Nartea    2 1 1 | 1 0 0  1 L2  L2  L1  2 1 1 | 1 0 0   2 3  b)  1 2 3 | 0 1 0    0 L3  L3  L1 5 7 1  |  2 1 0   L3  L3  L2   3 1 1 | 0 0 1   2 2 2     5 5 3  0  |  0 1  2 2 2   2 1 1 | 1 0 0  2 1 1 | 1 0 0   L  2 L   L  L  7 L 0 5  7 |  1 1 0   2 5 2  0 1  7 |  1 2 0   2 2 5 3   2 2 2   5 5 5  0 0 1 | 1 1 1  0 0 1 | 1 1 1     2 2  2 1 1 | 1 0 0 2 0 0 | 5 0 5      0 1 0 |  8 1 7  L1  L1  L2  L3  8 7  L1  12 L1   0 1 0 |  1    5 5  5 5  0 0 1 | 1 1     1  0 0 1 | 1 1 1     1 1  1 1 1 0 0 | 5 0 5   5 0 5      0 1 0 |  8 1 7 Deci A    1 8 1 7  5 5   5 5      0 0 1 | 1 1 1   1 1 1       Calculul rangului unei matrice prin metoda transformarilor elementare Se efectueaza transformari elementare asupra matricei pana cand toate elementele devin nule cu exceptia unor elemente de pe diagonala principala care devin unu Rangul matricei este numarul elementelor 1 de pe diagonala principala Aplicatie Determinati rangul matricelor 1 2 3 1 2  3 12    a) A    ; b) A    ; c) A   2 1 1 3 4 2 8   1 3 2    Solutie  1  1 2  L3 L3 3L1  1 2  L2   2 L2  1 2  C2 C2 2C1  1 0  a)             rang A=2 3 4  0 2  0 1 0 1 Lector univ dr Cristina Nartea  3 12  L2 L2  23 L1  3 12  L1 13L1  1 4  C2 C2 4C1  1 0  b)             rang A=1 2 8  0 0  0 0 0 0  1 2 3  L L 2 L  1 2 3  5   L32 L32  L1 1   L3  L3  3L2 c)  2 1 1    0 3  7     1 3 2  0 5 5        L2   1 L2 1 1 2 3 2 3   3   L3   20 L3   3   7  C2 C2  2C1 0 3 7    0 1     3 20  0 0 0    0 1   3  1 0 3 1 0 0     7 1 0 0 0 7  C3 C3 3C1  7  C3 C3  3 C2   1  0 1    0 1 0   rang A =3  3  3 0 0 1 0 0 1 0 0 1      in practica, pentru rezolvarea unui sistem de ecuatii liniare, procedam astfel: se efectueaza transformari elementare asupra matricei extinse pana cand toate elementele de sub diagonala principala devin nule Pe parcursul algoritmului pot aparea urmatoarele situatii:  coeficientii unei ecuatii devin toti nuli, iar termenul liber corespunzator este nenul, caz in care sistemul este incompatibil;  coeficientii unei ecuatii sunt toti nuli si termenul liber corespunzator este nul, atunci ecuatia respectiva este consecinta a celorlalte (deci inutila) Metoda transformarilor elementare consta in reducerea sistemului (1 2) la un sistem mai simplu, urmand pasii Pasul 1 Se schimba ecuatiile intre ele astfel incat prima necunoscuta x1 are coeficientul nenul in prima ecuatie, adica a11  0 Pasul 2 Pentru fiecare i>1, se aplica operatia Li  ai1L1  a11Li Lector univ dr Cristina Nartea Adica se inlocuieste ecuatia i cu ecuatia obtinuta din inmultirea primei ecuatii cu ai1 , inmultirea celei de a i ecuatii cu a11 si adunarea acestora Se obtine asfel o forma echivalenta a sistemului (i e are aceeasi solutie) a11 x1  a12 x2    a1n xn  b1   a2' j2 x j 2    a2' n xn  b2'     ' amj x j 2    amn ' xn  bm'  2 x j2 este prima necunoscuta cu coeficient nenul dintr-o alta ecuatie in afara de prima Se continua procedeul pana cand se ajunge la forma echivalenta a11 x1  a12 x2    a1n xn  b1   a2 j2 x j 2    a2 n xn  b2  (1 3)    arjr x j r    arn xn  br  Am notat coeficientii cu aceleasi litere ca in sistemul (1 2), dar in mod evident ei reprezinta alti scalari x  2 y  3 Aplicatia 1 Sa se rezolve sistemul  3x  6 y  9 Solutie  1 2 | 3  L2 3 L1  1 2 | 3       3 6 | 9 0 0 | 0 Deoarece ultima linie este numai cu zerouri, inseamna ca ecuatia corespunzatoare este inutila in acest caz, r=1<2=n, deci sistemul este compatibil nedeterminat cu 2-1=1 necunoscute secundare Fie aceasta y Atunci solutia sistemului are forma 3  2 y, y  , y   Lector univ dr Cristina Nartea  x  2 y  4 Aplicatia 2 Sa se rezolve sistemul  3x  6 y  9 Solutie  1 2 | 4  L2 3 L1  1 2 | 3        3 6 | 9   0 0 | 3  Sistemul este incompatibil 2 x  y  2 z  3w  1  Aplicatia 3 Sa se rezolve sistemul 3x  2 y  z  2w  4 3x  3 y  3z  3w  5  Solutie  2 1 2 3 | 1  L 3 L  2 L  2 1 2 3 | 1    L32 3 L11  2 L32    3 2 1 2 | 4     0 1 4 5 | 5   3 3 3 3 | 5   0 3 12 15 | 7       2 1 2 3 | 1   L3 3 L2  L3     0 1 4 5 | 5  Sistemul este deci incompatibil  0 0 0 0 | 8     x  2 y  3z  4  x  3 y  z  11  Aplicatia 4 Sa se rezolve sistemul  2 x  5 y  4 z  13 2 x  6 y  2 z  22 Solutie 1 2 3 | 4 1 2 3 | 4 1 2 3 | 4   L1  L2   LL32 2 L1  L3   L3  L2  L3   1 3 1 | 11  L4 2 L1  L4  0   1 4 | 7  L4 2 L2  L4  0   1 4 | 7 2 5 4 | 13  0 1 2 | 5 0 0 2 | 2       2 6 2 | 22  0 2 8 | 14  0 0 0 | 0 Lector univ dr Cristina Nartea  x  2 y  3z  4  Deci, sistemul este echivalent cu  y  4 z  7 Din ultima ecuatie, aflam z=1 inlocuind in a  2z  2  doua ecuatie, se obtine y=3, iar apoi, din prima ecuatie x=1 Deci, sistemul are solutie unica, iar solutia sistemului este 1,3,1  x  2 y  2 z  3w  2  Aplicatia 5 Sa se rezolve sistemul 2 x  4 y  3z  4w  5 5 x  10 y  8 z  11w  12  Solutie  1 2 2 3 | 2  L 2 L  L  1 2 2 3 | 2    L32 5 L11 L32    2 4 3 4 | 5     0 0 1 2 | 1   5 10 8 11 | 12   0 0 2 4 | 2       1 2 2 3 | 2   L3 2 L2  L3     0 0 1 2 | 1  0 0 0 0 | 0   Deci r=2<4=n Avem 4-2=2 necunoscute secundare Sistemul are forma echivalenta  x  2 y  2 z  3w  2   z  2w  1 Solutia sistemului este  4  2 y  w, y,1  2w, w , y, w   interpretare geometrica Pentru un sistem de doua ecuatii cu doua necunoscute a1 x  b1 y  c1  a2 x  b2 y  c2 pot aparea trei situatii a) Sistemul este incompatibil; b) Sistemul are solutie unica; Lector univ dr Cristina Nartea c) Sistemul este compatibil nedeterminat Reprezentarea grafica in plan a unei ecuatii de forma ax  by  c este o dreapta interpretarea geometrica a situatiilor de mai sus este a) Cele doua drepte sunt paralele; b) Cele doua drepte se intersecteaza intr-un singur punct; c) Cele doua drepte coincid Tema Sa se rezolve sistemele  x  2 y  3z  1  a) 3x  y  2 z  7 ; 5 x  3 y  4 z  2  2 x  y  2 z  10  b) 3x  2 y  2 z  1 ; 5 x  4 y  3z  4   x  2 y  3z  6  c) 2 x  y  4 z  2 4 x  3 y  2 z  14   x  3 y  4 z  2w  5  d)  2 y  5z  w  2  y  3z 4  